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In the limit of infinite spatial dimensions a thermodynamically consistent the- 
ory, which is valid for arbitrary value of the Coulombic interaction {U < oo), 
is built for the Hubbard model when the total auxiliary single-site problem ex- 
actly splits into four subspaces with different "vacuum states" . Some analytical 
results are given for the Hartree-Fock approximation when the 4-pole structure 
for Green's function is obtained: two poles describe contribution from the Fermi 
liquid component, which is ferromagnetic and dominant for small electron and 
hole concentrations ("overdoped case" of high-Tc's), whereas other two de- 
scribe contribution from the non-Fermi liquid, which is antiferromagnetic and 
dominant close to half filling ("underdoped case"). 
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1. Introduction 

In the last decade the essential achievements in the theory of the strongly corre- 
lated electron systems are connected with the development of the dynamical mean- 
field theory (DMFT) proposed by Metzner and Vollhardt [Q for the Hubbard model 
(see also Refs. [0, Q and references therein). There are no restrictions on the \J value 
within this theory and it turns out to be useful for intermediate coupling {U ~ t) for 
which it ensures the correct description of the metal-insulator phase transition and 
determines the region of the Fermi-liquid behavior of the electron subsystem. More- 
over, some classes of binary-alloy-type models (e.g., the Falicov-Kimball model) can 
be studied almost analytically within DMFT But in the case of the Hubbard 
model, the treatment of the effective single impurity Anderson model is very com- 
plicated and mainly computer simulations are used, which calls for the development 
of the analytical approaches 0. 
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Such approaches can be built using a systematic perturbation expansion in terms 
of the electron hopping ^ using a diagrammatic technique for Hubbard operators 
0, H]. One of them was proposed for the Hubbard {U = oo limit) and t — J models 
. The lack of such approach is connected with the concept of a "hierarchy" system 
for the Hubbard operators when the form of the diagrammatic series and the final 
results strongly depend on the system of the pairing priority for Hubbard operators. 
On the other hand it is difficult to generalize it to the case of arbitrary U. 



In our previous paper ||T0[ we developed for the Hubbard-type models a rigorous 
perturbation theory scheme in terms of electron hopping that is based on the Wick 
theorem for Hubbard operators is valid for arbitrary values of f/ (f/ < oo) 

and does not depend on the "hierarchy" system for X operators. In the limit of 
infinite spatial dimensions, this analytical scheme allows us to build a self-consistent 
Kadanoff-Baym type theory for the Hubbard model and some analytical re- 
sults are given for simple approximations. Here we shall consider possible magnetic 
orderings in the Hartree-Fock type approximation. 

2. Perturbation theory in terms of electron hopping 

We consider the lattice electronic system described by the statistical operator: 

p = e-^^»a(/3), a(/3)=Texpi-|rfr|rfr'5^tr^.(r-r')aL(r)a,.(r')L (2.1) 

I J 

where Hq = Yli Hi is a sum of the single-site contributions and for the Hubbard 
model we must put 

Hi = Uni^riii - fi{ni^ + riii) - hlui-i - riii), ^iji'^ - '^') = UjS{T - t'). (2.2) 

It is supposed that we know the eigenvalues and eigenstates of the zero-order 
Hamiltonian: Hi\i,p) = \p\i,p), and one can introduce Hubbard operators Xf^ = 
\i,p){i,q\ in terms of which the zero-order Hamiltonian is diagonal 

Ho = J2J2^P^f"- (2-3) 

i p 

For the Hubbard model we have four states \i,p) = |2,nj|,nj|): \i,0) = \i,0,0) 
(empty site), \i,2) = |i, 1, 1) (double occupied site), = \i,l,0) and \i, I) = 

\i, 0, 1) (sites with spin-up and spin-down electrons) with energies Aq = 0, A2 = 
U — 2/i, A| = h — fi, and A^ = —h — fi. The connection between the electron 
operators and the Hubbard operators is the following: 

n,, = Xf + Xr; a,, = X^," + aXp. (2.4) 

The expression for {a{l3))o is a series of terms that are products of the hopping 
integrals and averages of the electron creation and annihilation operators or Hubbard 
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operators, that are calculated with the use of the corresponding Wick's theorem 
0, 1^ , and can be written as fl^ : 



exp 





where arrows denote the zero-order Green's functions 

1 



9pq{ 



(2.5) 





(2.6) 



wavy lines denote hopping integrals and □, . . . stand for complicated "n vertices" . 
Each vertex (many-particle single-site Green's function) is multiplied by a diagonal 
Hubbard operator denoted by a circle. 



3. Irreducible many-particle Green's functions 

For the Hubbard model, expressions for the two-vertex 

® =5^XfV(p)(^.), (3.1) 
p 

the four-vertex 

] [ = X^^g„(p) {uJn ) ga(p) {uJn+m.)Uaa (p) {^n,^l\^m.) 9a{p) {^l ) 9a{p) j^l+m), (3.2) 



and for the vertices of higher order possess one significant feature |TD[. They de- 
compose into four terms with different diagonal Hubbard operators X^^, which 
project our single-site problem on certain "vacuum" states (subspaces), and zero- 
order Green's functions 

which describe all possible excitations and scattering processes around these "vac- 
uum" states. Here 

n , , I, , ^ - / ^ ± U^920i^n+l+ra) ioi p=0,2 

f/.,(p)(^„, uilum) - I ^ ^ U'g^,iuj^^_i) for p = a,a ^^'^^ 

is a renormalized Coulombic interaction in the subspaces. In diagrammatic notations 
expression (|3^ ) can be represented as 

n - X ^ >-< . (3.5) 
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where dots denote the Coulombic correlation energy U and the dashed arrows denote 
bosonic zero-order Green's functions: doublon (72o(^m) or magnon g^^^Um)- The 
contributions to the six-vertex can be presented by the following diagrams: 



^ -^-r- ^rr^ >-^-< 



(3.6) 



where the first three diagrams contain the internal vertices of the same type as in 
Eq. ( |3.5| ). So, we can introduce primitive vertices by which one can 

construct all n vertices in the expansion ( p.5|) .H 

4. Dynamical mean-field theory 

Within the framework of the considered perturbation theory in terms of electron 
hopping a single-electron Green's function can be presented in the form 

where we introduce an irreducible part Ho-(i-o'„, k) of Green's function which, in gen- 
eral, is not local. In the case of infinite dimensions d oo one should scale the 
hopping integral according to t.y tij/ \fd in order to obtain finite density-of-states 
and it was shown by Metzner in his pioneer work [|l^ that in this limit the irreducible 
part become local "Eijaij — t') = 6ijE(j{T — r') or Ho-(co'„,fc) = E^j^Un) and such a 
site-diagonal function can be calculated by mapping the d ^ oo lattice problem 
( |2.1| ) onto the atomic model with the auxiliary Kadanoff-Baym field 



The self-consistent set of equations for So-(ti;„) and Jo-(co'„) (e.g., see Ref. and 
references therein) is the following: 

^ E \ , = / w ^ = ^i'^i^n, {J.M}), (4.3) 

where G^^\uJn, {Jai^^n)}) is the Green's function for the atomic hmit ( [4.2| ). 

The grand canonical potential for the lattice is connected with the one for the 
atomic limit by the expression 0] 

% = ^a-iy:< InG^'^n^n) - ^ > :inG.(^„, k) } . (4.4) 



^ riCT I k J 



On the other hand, we can write for the grand canonical potential for the atomic 
limit Qa the same expansion as in Eq. ( p.5|) but with diagonal X operators at the 

^For n vertices of higher order a new primitive vertices can appear but we do not check this 
due to the rapid increase of the algebraic calculations with the increase of n. 
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same site. We can reduce their product to a single X operator that can be taken 



outside of the exponent in (|2.5| ) and its average is equal to (Xpp)o = e e ^^"^ 

Finally, for the grand canonical potential for the atomic limit we get [III 



^ln^e-^^(-), (4.5) 



where ^(p) are the "grand canonical potentials" for the subspaces. 

Now we can find the single-electron Green's function for the atomic limit by 

G'^:\r - r') = jj^^^ = E - r'), (4.6) 

where G'cr(p) {t — t') are single-electron Green functions for the subspaces character- 
ized by the "statistical weig hts" Wp = e-^^w /X;, e-^^(9) and our single-site atomic 
problem exactly splits into four subspaces p = 0,2, |, |. 

The fermionic Green's function in subspaces can be written as 

Ga{p)iuJn) = ^ . ^_ J f — y -a{p)M = + /i<x " ^^^(p) " S^(p)(u;„), (4.7) 



where 'T-^lp) = —dXp/d^^ = for p = 0, a and 1 for p = 2, o" is an occupation of the 
state \p) by the electron with spin a, and the self-energy So-(p)(ci;„) depends on the 
hopping integral J^iiun') only through quantities 

'^a'(p){^n') = Ga'{p){uJn') — '^a' {p){^n') ■ (4.8) 

Now, one can reconstruct expressions for the grand canonical potentials in 
subspaces from the known structure of Green's functions: 



^In (1 - J^{uJn)Ea{p)M) -^Yl Sa(p)(Wn)^a(p)(^n) (4.9) 



where $(p) is some functional, such that its functional derivative with respect to 
^ produces the self-energy: (5$(p)/5\E'o-(p)('^n) = ^a{p){^n)- So, if one can find or 
construct the self-energy So-(p) i^n) he can find Green's functions and grand canonical 



potentials for the subspaces and, according to Eqs. ( [4.5|) and (|4.6| ), can solve atomic 
problems. 

From the grand canonical potential ( [4.5| ) and ( [4.9| ) we get the following mean 
values 

na = J2 V-(P)^ ^<x{p) = ^f[p) + ^<x(p)(^n) - (4.10) 



where in the last term the partial derivative is taken over not in the chains ( [4.8|) . 

For the Falicov-Kimball model JiioJn) = and, as a result, S|(p)(u;n) = and 
S|(p)(co'„) = g^(^j,){(jjn) which immediately gives results of Ref. @] (see also Ref. 

For the Hubbard model there is no exact expression for the self-energy but the 
set of Eqs. (|4.7|) and (|4.9| ) allows one to construct different self-consistent approxi- 
mations. 
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5. Hartree-Fock approximation 

One of the possible approximations is to construct the equation for the self-energy 
in the following form: 

Sa(p)K) = -n ^U'^a{p){uJn'), (5.1) 
^ n' 

which, together with the expression for mean values 

,(0) 



gives for the Green's functions in the subspaces an expression in the Hartree-Fock 
approximation: 



Ga{p){u}n) = — — — 77Z jr. \ - (5-3) 



1 

Now, the grand canonical potentials in the subspaces are equal 



fi{p) = Ap- X In [1 - J,(a;„)S,(p)(^„)] - U (n^^p^ - (n^^p) - nf^^^ (5.4) 

ncr 

and the Green's function for the atomic problem (|4.6| ) yield a four-pole structure, 
that, in contrast to the alloy-analogy solution So-(p)(a;„) = 0, possesses the correct 
Hartree-Fock limit for small Coulombic interaction U <^ t. On the other hand, in 
the same way as an alloy-analogy solution, it describes the metal-insulator transition 
with the change of U . 

In Ref. [0] it was shown that the main contributions into the total spectral 



weight function Pa{uj) = ^'^Ga {uj — iO"*") come from the subspaces p = for the 
low electron concentrations (n < |, /i < 0), p = 2 for the low hole concentrations 
(2 — n < |, yU > f/) with p = a,a for the intermediate values. For the small electron 



or hole concentrations, the Green's function for the atomic problem ( [4.6| ) possess the 
correct Hartree-Fock limits too. It was supposed that the Hubbard model describes 
strongly-correlated electronic systems that contain four components (subspaces). 
Subspaces p = and p = 2 describe the Fermi-liquid component (electron and hole, 
respectively) which is dominant for the small electron and hole concentrations, when 
the chemical potential is close to the bottom of the lower band and top of the upper 
one ("overdoped regime" of high-Tc's). On the other hand, subspaces p =t and | 
describe the non- Fermi-liquid (strongly correlated, e.g., RVB) component, which is 
dominant close to half-filling ("underdoped regime"). 

Here we consider possible magnetic orderings. At low temperatures the p = and 
p = 2 components for the low electron or hole concentrations are in the ferromagnetic 
state, while the non-Fermi-liquid one is antiferromagnetic (AF) close to half-filling, 
see Fig. |l]. For the intermediate concentration values the picture is very complicated, 
even frustrated. It is due to the fact that equations for the mean values ( |5.2| ) have 
several solutions in this region, which, on the other hand, are mutually connected 
with the dynamical mean field Ja{oJn)- It is difficult to determine the ground state 
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n 2.0 



Figure 1. Ferromagnetic mp and antiferromagnetic niAF order parameters vs 
electron concentration for U = 1.56, T = 0.14. 



for this, possibly "pseudo-gap" , region, which is located between the ferromagnetic 
and antiferromagnetic phases. 

In Fig. ^ we presented the phase diagram (T, U) — the temperature of the AF 
ordering vs correlation energy U, which is in a qualitative agreement with the results 
of Refs. and reproduces the results of the Hartree-Fock theory and mean 
field approximation for U ^ t and U ^ respectively. Our results for the AF 
critical temperature for small U are higher then the one of the Quantum Monte 
Carlo simulations fl^ by about a factor at three that describes the reduction of the 



Hartree-Fock solution by the lowest order quantum fluctuations |16 




Figure 2. Phase diagram (T, U) at half- filling n = 1 (AF — antiferromagnetic 
phase, PM — paramagnetic phase). 
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